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I. INTRODUCTION 

Rotating spiral waves (RSWs) are typical robust self- 
organized structures (autowaves) in various dissipative 
systems, ranging from chemistry, biology and ecology 
fig . to nonlinear optics and lasers [i^ • 

Evolution of these spatio-temporal structures is now 
a subject of extensive theoretical studies, experimental 
investi gati ons and computer modeling (see, e.g., [TqIIs^ 
l4(i lil . I43 and references therein) . 

Many mathematical models in this field are based on 
concept of excitability 35i|. An important approach in 
modern computer investigations of excitable systems is 
based on using of discrete parallel models with local in- 
teractions, mostly by cellular automata (CA). 

In the framework of the CA approach 0,1131 j the three- 
level representation of the excitable media states is an 
effective tool for carrying computer experiments with 
multi-particle systems possessing long-time nonstation- 
ary dynamics (see 28.] and references therein). Paral- 
lelization of updating of the active medium states and 
locality of interactions between elementary parts of the 
medium are the fundamental principles of the CA ap- 
proach. Using of this approach admits to fulfill modeling 
of large systems by emulation of parallel matrix transfor- 
mations. 
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The essence of three-level representation of excitable 
medium |^ [s^ is as follows. Each elementary part 
or active center (AC) of an excitable system has the single 
stable ground state, say Li, and at least two metastable 
(upper) states, say Lu and Lni- 

The higher metastable state Lm is the excited one, 
it may be reached only after certain hard perturbation 
(from neighboring AGs or from external source) of the 
ground state L\ of the considered AC. A group of excited 
ACs may excite another ACs etc. Due to metastabilty, 
the lifetime Te of the excited state Lm is a finite value, 
after which every excited AC spontaneously reaches the 
intermediate refractory state Lu. 

Refractorily means "sleeping" state at which no excita- 
tion of the considered AC is possible both from neighbor- 
ing ACs or from an external source. Moreover, the AC in 
refractory state can't participate in excitation of another 
ACs (lying at Li). So, after a refractorily lifetime Tr, the 
considered AC reaches the ground stale Lj. Usually only 
the cyclic transitions (Lj — > Lm Lu Li ...) are 
permitted. 

The concept of excitability was developed in biology, 
but it is widely used now in chemistry and physics. In 
particular, this concept have been applied to solve many 
problems arisin g in modern nonlinear optics, laser physics 

etc. [13 ini [Till El El 113 

In '2^ we have considered even more close laser-like 
analog of excitable system, namely the three-level ac- 
tive medium of the microwave phonon laser (phaser) with 
dipole-dipole interactions between AC. The only signifi- 
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cant difference of phaser medium from the usual excitable 
one is the second channel of diffusion of excitations. 
This system was experimentally studied in Institute 
of Radio-Physics and Electronics (Kharkov, Ukraine) 
[22, 23^ 2^ i2hi 2^ i22ji demonstrating self-organization, 
bottlenecked cooperative transient processes and other 
nonlinear phenomena under extremely low level of intrin- 
sic quantum noise (phaser has 15 orders lower intensity 
of spontaneous emission in comparison to usual optical- 
range lasers). 

In preceding publications [2^ |2^ [s^, |^ we carried 
out a series of computer experiments on RSWs dynam- 
ics in autonomous excitable systems with various control 
parameters and random initial excitations. The most in- 
teresting phenomena of self-organized vorticity observed 
in computer experiments .28„ .29, ,,30». <3Lj were as follows: 
(a) spatio-temporal transient chaos in form of highly 
bottlenecked collective evolution of excitations by RSWs 
with variable topological charges; (b) competition of left- 
handed and right-handed RSWs with unexpected fea- 
tures, including self-induced alteration of integral effec- 
tive topological charge; (c) transient chimera states, i. e. 
coexistence of regular and chaotic domains in excitable 
media; (d) branching of an excitable medium states with 
different symmetry which may lead to full restoring of 
symmetry of imperfect starting pattern. Phenomena (a) 
and (c) are directly related to microwave phonon laser 
dynamics features observed earlier in real experiments at 
liquid helium temperatures on corundum crystals doped 
by iron-group ions 22, 23, 24, 25, 26, 27J. 

In the present paper, we report some qualitatively 
new features of coexistence, competition and collapse of 
RSWs in three-level model 28, 29, 30, 31] of excitable 
media under conditions of strong influence of the sec- 
ond channel of diffusion of excitations. Part of these 
features are caused by unusual mechanism of RSWs evo- 
lution when RSW's cores get into the surface layer of 
an active medium (i. e. the layer of AGs resided at the 
absorbing boundary). Instead of well known scenario of 
RSW collapse, which takes place after collision of RSW's 
core with absorbing boundary, we observed complicated 
transformations of the core leading to "reflection" of the 
RSW from the boundary or even to birth of several new 
RSWs in the surface layer. To our knowledge, such non- 
linear "reflections" of RSWs and resulting die hard vor- 
ticity in excitable media with absorbing boundaries were 
unknown earlier. 



II. SOME REMARKS ON ZYKOV-MIKHAILOV 
MODEL AND ITS TWO-CHANNEL 
MODIFICATION 

The Zykov-Mikhailov (ZM) model [111 is a discrete 
parallel mapping with local interactions, which may be 
defined as three-level CA Let an active discretized 
medium Tic has the form of rectangular 2D-lattice. Each 
cell of the lattice contains the single three-level AC with 



coordinates (i,j), where (min(i) = 1), (min(j') = 1). All 
the AGs in the Tie are identical, and they interact by the 
same set of rules (ZM model, as well as most of the CA 
models, is homogeneous and isotropic in the von Neu- 
mann sense fsSj)- Each AG has the single stable ground 
state Li, and two metastable states Lu and Lm. 

The upgrade of states S^f = ^("Hi, j) of aU AGs (i. e. 
their intralevel or interlevel transitions) are carrying out 
synchronously at each step < n < during the system 
evolution. All intralevel (Lj — > Li, Lu — > Lu, Lni — > 
Lni) and some of interlevel (namely Li — > Lni, Lui — > 
Lu, Lu Li) transitions are permitted. The final step 
N may be either predefined or it will be searched during 
CA evolution as time (i.e. quantity of discrete steps) 
for reaching an attractor of this CA. The conditions of 
upgrade [^H^ depend both on sl^~^^ and on S'|,"7"^\ 
where belongs to certain active neigborhood of 

AG at site In the framework of the ZM model 

(and in most of other CA models of excitable chemical 
systems) , diffusion of excitations is possible only by Lj — > 
Lui transition. This is a single-channel (IC) mechanism 
of diffusion. 

A modification of ZM model was proposed by us 
[28L l2^ I30I l3lj| to adapt it for emulation of some as- 
pects of dynamics of class-B optical lasersJE^JS [s^l and 
microwave phonon lasers [23, 123, l24L l25l l2a. l27j . The re- 

laxational properties of AGs in microwave phonon lasers 
are of the same type, as in class-B optical lasers fisi . 
and there are experimental confirmations [23L 123 . 125ll26|| 
of common properties of microwave phonon lasers and 
class-B optical lasers. 

In order to formulate boundary conditions for a Tie of 
finite size (i. e. when max(j) = Mj5c,max(j) = My), a 
set Tlv of virtual cells with coordinates (i 0) V {j ~ 
0) V a = Mx + 1)V {j = My + 1) may be introduced 
|2Ml29ll30ll3ll |. Each virtual cell contains the unexcitable 
center having the single level Lq. The surface layer of a 
bounded active system is represented by the set OTsurf of 
AGs with coordinates (i = 1) V {j = 1) V (j = Mx) V 
(j = My). Note that our surface layer is of ID type 
because active medium itself is of 2D type (we stress this 
circumstance to avoid misunderstandings). 

The original ZM model has, as it was pointed out 
earlier, the single channel of diffusion |5^. Such 
IG-diffusion models are adequate for chemical reaction- 
diffusion systems [4^ . In microwave phonon laser active 
system, the multichannel diffusion of spin excitations is 
the typical case, because it proceeds via (near)-resonant 
dipole-dipole (d — d) magnetic interactions between para- 
magnetic ions ,22]. For a three- level system, which is 
the simplest microwave phonon laser system, there are 
3 possible channels of resonant diffusion. In the case 
when d — d interactions are forbidden at one of three 
resonant frequencies of a three-level system, the asym- 
metric two-channel (2G) diffusion is realized under con- 
ditions of perfect refractority of the intermediate level — 
see jl^ I23, 113, 1^. Note that asymmetric diffusion is 
well known in biology (see, e.g., the book of D. A. Frank- 
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Kamenetzky [9j). Recently a kind of asymmetric diffu- 
sion was proposed by N. Packard and R. Shaw j3^J for a 
mechanical system. 



and 2,^"^ for 



Lpf^^ ) and two types of partial attributes uf^ ^^^^ 
each individual AC in Tic- Full description of all AG's 
possible states is as follows: 



III. THREE-LEVEL MODEL OF EXCITABLE 
SYSTEM WITH TWO-CHANNEL DIFFUSION 
(A MODIFIED ZYKOV-MIKHAILOV MODEL) 

A. States of active centers and branches of 
evolution operator 




if n = 0, 
if n ^ 0; 



(1) 



(2) 



Let dJle is a rectangular 2D lattice containing Mx x My 
three-level AGs. The upgrade of state 5*1"^ of each AG is 
carried out synchronously during the system evolution. 
The excited (upper) level Lni has the time of relaxation 
Te > 1, and the refractory (intermediate) level Ln has 
the time of relaxation > 1. Both Tg and are inte- 
ger numbers. This model of excitable system is a kind 
of three-level GA. The original ZM model [H Ell was 
formulated as three- level GA with IG diffusion. 

In our three-level model of excitable system (TLM), 
the second channel of diffusion of excitations is included. 
A detailed argumentation for using of the 2G diffusion 
mechanism one can find in , j8 j . 

In the TLM, the first channel of diffusion accelerates 
the transitions Li Lm for a given AG, and the second 
channel of diffusion (which is absent in the ZM model) 
accelerates the transitions Lm Lu. The complete 
description of AG's states 5*1"' in such the TLM model 
includes one type of global attributes (the phase counters 



S^\Lu) = {v^'). (3) 

In the framework of the TLM model, the phase coun- 
ters lie in the interval ip^J^ € [0, Tg + r^]. The following 

correspondences between cp^^"^ and Lk take place (by def- 
inition) for all the AGs in 071c at all steps n of evolution: 

{LK=U)^{^\f^Oy, (4) 
(LK-iiii)^(0<(/9|;^ <re); (5) 
{Lk = Lu) ^ (re < ip^f < t, + r,) . (6) 

These correspondences are of key significance for all 
multistep-relaxation models (Reshodko model 0, 0, j 
ZM model [^E^ etc.): there are only three discrete lev- 
els, and relaxation of each AG is considered as intralevel 
transitions (or transitions between virtual sublevels). 



The evolution of each individual AG proceeds by sequential cyclic transitions Lk Lk' (where K and K' G 
{I, III, II}), induced by the Kolmogorov evolution operator the TLM model, the evolution operator VL has 

three orthogonal branches fJi, Sim and f2ii, which we call ground, excited and refractory branches respectively. The 
choice of the branch at iteration n -I- 1 is dictated only by the global attribute of the AG at step n, namely: 

pi (4")), if4") = o, 

= n (^If) = nm ) , if < < T,, (7) 



B. Ground branch of the evolution operator 

At step 71 + 1, the branch Slj by definition fulfills operations only over those AGs, which have Lk = Li at step n. 
These operations are precisely the same, as in ZM model 0,11^, namely: 




p,q 
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where A^j'' is the accumulating term for the w-agent, which is an analog of chemical activator |2 



D<if is the 



first-channel diffusion term; h is the threshold for the w-agent {h > 0); g is the accumulation factor for the u-agent 
{g £ [0, 1]); C{p, q) is the active neighborhood of the AC at site («, j); and the u-agent arrives to AGs with Lk = Li 
only from AGs with Lk = Lm in C{p, q): 



t(") 



0, if(^(") 



j+q 
> r< 



(10) 



The definition of the diffusion term Dl^ in Eqn. © is very flexible. Apart of well-known neighborhoods of the 
Moore {C{p, q) = Cm{p, <?)) and the von Neumann (C(p, q) — Cn{p, q)) types: 



Cm{p, q) 



1, if 



( bl < 1) A ( \q\ < 1) A (<5po<550 ^ 1) 



0, otherwise, 



(11) 



1, if [(H = l)©(|,s| = l)J; 
0, otherwise, 



(12) 



one can easy define another neighborhoods, e. g. of the box type Cb{p, q) or of the diamond type Cd{p, q) 8] (which 
are the straightforward generalization of the Moore Cm{p, q) and the von Neumann Cn{p, q) neighborhoods), etc. 

On the other hand, the ZM definition 0,15^ of the weight factors Jj^"^ -j^^ HlOfl may be extended out of the binary 
set {0, 1} to model various distance-dependent phenomena within, e. g., Cd(jp, q), Csip, q)- 

In this work, however, we restricted ourselves by the Moore neighborhood and by weight factors of the form 
(|10(l . Another types of neighborhoods, weight factors and some other modifications of the model will be studied in 
subsequent papers. 



C. Excited branch of the evolution operator 



At the same step n + 1, the branch Om fulfills operations only over those AGs, which have Lk = Lm at step n: 



n+l 




if 



(0<(^, 



< r.) A < /) V(^, 

if (0 < ^ < T,) A 



(n) 
i'j 



> 



(13) 



(n+l) 



r^(^S.+,)-E^(^^'9)Q!:u,. (14) 



P,9 



where D.^^ is the second- channel diffusion term; / is the threshold for the z-agent (/ > 0), and we assume that 
z-agent arrives to excited AG {Lk = Lm) only from those AGs, which have Lk = Li in C{p, q): 

One can see from (|14ll that the z-agent does not accumulate during successive iterations. In other words, the branch 
fiiii at step n -I- 1 produces "memoryless" values of partial attributes ^i"'*'^'' for AGs having Lk = Lm at step n (in 
contrary to the w-agent for AGs having Lk — Li at step n). The z-agent may accelerate transitions from excited AGs 
to refractory ones. This is the most important difference between our 2G model of three-level excitable medium and 
the original IG model of ZM [Hi^. 
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D. Refractory branch of the evolution operator 

The branch flu does not produce/change any partial attributes at aU (because the intermediate level Lu is in the 
state of refractority) . It fulfils the operations only over those AGs, which have = Lu at step n: 

^^"^'^ (16) 

Generally speaking, there are many examples of active media with weak refractority (when the unit at the inter- 
mediate level Lu is not absolutely isolated from its neighbor hooding units). But in this work we restrict ourselves by 
the case of perfect refractority lit) I) , which is valid, e. g., for the microwave phonon laser systems of the Ni^+ : AI2O3 
type 




E. Reduction of the TLM model to the ZM automaton 



The second-channel diffusion gives the contribution to the TLM dynamics if / < A/", i.e. if / < 4 for C = Cat, 
/ < 8 for C = C]\f and so on. 

At ((C = Cat) a (/ > 4)) our TLM model is of ZM-hke (i.e. IG) type, and at ((C = Cm) A (/ > 8)) it becomes 
equivalent to the original ZM model ,21jj5^], which (with slight rearrangement of cases) is as follows: 



n+l 



(") I 1 



if < (^(;) < T, 



o)aO 



,("+!) 



(^1;) = 0) A «+^' <h) y{ 



(n+l) 



(17) 



where is defined as in lPl- l(TT)l . 



F. Geometry of active media, boundary conditions 
and transients in TLM 



Pattern evolution is very sensitive to geometry of active 
medium and boundary conditions even in the framework 
of operation of the same GA. There are three types of 
geometry and boundary condition (GBG) combinations 
most commonly used in computer experiments with GA 
(see Table I). Let us consider them in details. 

• The GBG-1 type is defined as follows: 



where 



Mx if iv = 0; 
1, if = Mx + 1; 
iv, otherwise, 



(18) 



(19) 



My ifjv=0; 
1, ifjv^My + l; 
jv, otherwise. 



(20) 



Here {iv,jv) G and the coordinates (io , j^) belong to 
the border set Tl^ C VJle of excitable area, i. e. (ic, jc) = 
if (i = 1) V (j = 1) V (i = Mx) V (j = My). 



• The GBG-2 type can be easily defined by using an 
extended interval for the phase counter (p. Let every cell 
(*v,Jv) in 2)Iv contains one unexcitable unit. All these 
unexcitable units have frozen (p(")(iv,Jv) atallnG [0,N] 
without reference to states of {ic,jc)j ^^d 



(n) 



(«v, jv) 



X — const. 



(21) 



where, x > Te -I- or x < 0. These unexcitable units 
in 3Jlv "absorb" excitations at the border of 9Jlc in con- 
trary to the case of GBG-1, where excitations at the bor- 
der 9Jlo are reinjected in active medium. In many cases 
this difference leads to qualitatively different behaviour 
of the whole TLM. Note that by this way one may also 
define TLM with inhomogeneous active medium, where 
some of AGs are changed to unexcitable units, i. e. un- 
excitable units are placed not only in SWv, but in SJte 
too. This may be easily done by introducing an addi- 
tional orthogonal branch D,q into the evolution opera- 
tor n. The pointed additional branch is activated at 



(n) 



= X only and is simply the identity operator fig — /, 



I.e.: tpy, ■ = ^of 



in) 



flj"^ . And even more complex 



behaviour of such "impurity" units may be defined using 
analogous approach (i. e. by introducing additional or- 
thogonal branches in an evolution operator): there may 
be pacemakers [2lj or another special units embedded in 
active medium and described by phase counter 7^ x in 
extended areas {(p > Te + Tr) or {ip < 0). 
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TABLE I: Geometry of active media and boundary conditions 



Type 


Geometry 


Boundary conditions 


GBC-1 


Bounded, Toroidal 


Cyclic 


GBC-2 


Bounded, Flat 


Zero-flow 


GBC-3 


Unbounded, Flat 


(free space) 



• The GBC-3 type is, strictly speaking, a case of bor- 
derless system. But the starting pattern has, of course, 
bounded quantity of AGs with Lk 7^ Li, located in 
bounded part of active medium. Hence, during the whole 
evolution (1 < n < N) the front of growing excited area 
will meet the unexcited (but excitable!) AGs only. For 
the case GBG-3, there is neither absorption of excitations 
at boundaries (as for GBG-2), nor feedback by reinjection 
of excitations in the active system (as for GBG-1). From 
this point of view GA of GBC-3 type are "simpler" , than 
GA of GBG-1 and GBG-2 types. On the other hand, CA 
of GBC-3 type are potentially infinite discrete systems 
where true aperiodic (irregular, chaotic) motions are pos- 
sible, in contrast to finite discrete systems, possessing, of 
course, only periodic trajectories in certain phase space 
after ending of transient stage [2(il |. 

On the other hand, there is also very important dif- 
ference between GBC-2 and other two types of GBC. 
The system of GBC-2 type is the single from the three 
ones under consideration which interacts with the ex- 
ternal world dynamically (besides of relaxation). Sure 
enough, GBG-1 and GBG-3 are connected to this outer 
world only by relaxation channels (the Tg and are the 
measures of this connection). In contrary, a system of the 
GBC-2 type interacts with surroundings through the real 
boundary, which is in fact absent in toroidal finite-size 
active medium of GBG-1 type and it is absent by defini- 
tion for a system of GBG-3 type. Despite of elementary 
mode of such interaction (boundary simply "absorbs" 
the outside-directed fiow of excitations from («o, jf), a 
system of the GBC-2 type may demonstrate very spe- 
cial behaviour. The main attention in this work is de- 
voted just to TLM of GBC-2 type as the most realistic 
model of microwave phonon laser active system, where 
both the mechanisms of interaction (dynamical and re- 
laxational) of an active medium with the outer world are 
essential. Dissipation in a microwave phonon laser active 
medium (highly perfect single crystal at liquid helium 
temperatures) is caused by two main mechanisms: (a) 
dynamical, by coherent microwave phonon and photon 
emission directly through crystal boundaries and (b) re- 
laxational, by thermal phonon emission. There are, of 
course, many more or less important differences between 
the TLM model and the real microwave phonon laser 
systems 0, HE j but in any case autonomous mi- 
crowave phonon laser has only these two mechanisms of 
interaction with the outer world. 



G. Initial conditions 

Patterns in ZM are described in terms of levels (or "col- 
ors") of AGs, and initial conditions are formulated simply 
as matrix of AGs levels Lx- But states of AGs in such 
the automata as ZM or TLM are fully defined not only by 
levels Lk itself. There are global attributes which must 
be predefined before TLM evolution is started. Some of 
partial attributes must be predefined too. These points 
are essential for reproducing of the results of computer 
experiments with ZM, TLM, etc. 

In our work such the initial conditions for the global 
attributes ipij are used: 

^^l, ifi^^-iiii; (22) 

Initial conditions for Uy must be defined for ground-state 
AG's [l'^k ^ -^i) oiily- III this work, we suppose 

(z.i?=0) IFF ^f^O), (23) 

where IFF means "if and only if" . Initial conditions for 
Zij are undefined for all because z-agent is not de- 
fined at n = 0, see Eqn. |(5J. As the result, initial con- 
ditions will be defined in our TLM model as the matrix 
(where (^l^-* G {0, l,Te + 1}) with additional con- 
dition given by Eq. (|23|l . 

IV. RESULTS AND DISCUSSION 

A bounded solitary domain of excited AGs having 
appropriate relaxation times and placed far from grid 
boundaries (or at unbounded grid) may evolve to RSWs 
if and only if there is an adjoined (but not a surround- 
ing) domain of refractory AGs 0, H, 0, • In this case 
RSWs apear by pair with opposite sgn{QT) and inte- 
gral topological charge = J2m <3t™^ obviously 
conserved (by infinity in time if grid is unbounded). If 
such solitary excited-refractory area is placed in a start- 
ing pattern near the grid boundary and the GBG-2 condi- 
tions take place, the single RSW appears and evolves and 
^(intcgi) .^^ course, not conserved in this case. These 
simplest and well known examples illustrate possibilities 
of coexistence and competition of one or two RSWs in 
excitable media, and possible scenarios of the RSW(s) 
evolution may be easily forecasted. 

An evolution of complex patterns with multiple, irreg- 
ularly appearing, chaotically-like drifting and colliding 
RSWs (or, possibly, another spatio-temporal structures) 
is in essence unpredictable without direct computing of 
the whole transient stage. So, the best way to investigate 
cellular automaton (TLM in particular) is to run it, be- 
cause, as S. Wolfram pointed out, "their own evolution 
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is effectively the most efBcient procedure for determining 
their future" (see page 737). 

Here we present results of our computer experiments 
with 2C model of excitable media described in Section 
mil The main tool used in these experiments was the 
cross-platform software package "Three-Level Model of 
Excitable System" (TLM) © 2006 S. D. Makovetskiy 
[3ll | . It is based on extended and improved algotithms of 
discrete modeling of three-level many-particle excitable 
systems, proposed by us in [29ll30ll3l|. T he TLM pack- 
age is written in the Java 2 language |l5l | using the Java 
2 SDK (Standard Edition, version 1.4.2) and the Swing 
Library. 

A discrete system with finite set of levels may have only 
two types of dynamically stable states at bounded lattice. 
The first of them is stationary state, and the second is 
periodic one. They may be called attractor s by an alog y 
with lumped dynamical systems (see e. g. fla. Isil ISq'). 
For our cellular automaton, the first of such the attrac- 
tors is spatially-uniform and time-independent state with 

^If = 0, where {i G [l,Mx]) A (j G [l,My]), n > nc- 
In other words, the only stationary state of TLM is the 
state of full collapse of excitations at some step nc (by 
definition of excitable system). The second type of TLM 
attractors includes many various periodically repeated 
states (RSW is a typical but not the single case) . In this 
case v?!;^ = where {i G [l,Mx]) A (j G [l,My]), 

n > np] np is the first step of motion at a periodic at- 
tractor; T is the integer-number period of this motion, 
T > 1. 

If starting patterns are generated with random spa- 
tial distribution by levels, then the time intervals nc or 
np may be considered as times of full ordering in the 
system. Such irregular transients are of special inter- 
est from the point of view of nonlinear dynamics of dis- 
tributed systems 1, 6, 14, 33, 34] because they may be 
bottlenecked by very slow, intermittent morphogenesis 
of spatial-temporal structures. In the next Subsections 
we study this collective relaxation for cases of collapse 
and periodic final states of TLM evolution, including an 
important case of lethargic transients. 

A. Collapse of RSWs 

In this Subsection, we describe a typical scenario of the 
two-dimensional, three-level and two-channel-diffusion 
(2D-3L-2C) excitable system evolution leading to the full 
collapse of excitations (when the system reaches the sin- 
gle point-like spatio-temporal attractor — stable steady 
state with Lk = Li for all AGs). Some typical stages of 
such the scenario of evolution are shown at Figure 1. 

The system evolves at first stages (0 < n < lO'') 
to a labyrinthic structure which contains several RSW- 
nucleating domains (Figure 1, n = 10000). 

During the next stage of the evolution, the large-scale 
vortex is formed (Figure 1, n = 60000). This vortex 
pushes the labirinthic structure to the boundaries of the 



active medium. As the result, the vortex occupies the 
whole active medium (Figure 1, n = 100000). 

But this giant rotating structure is unstable. Its core 
splits into RSWs with Qt = —1, and at this stage the 
active medium is filling by such the RSWs (Figure 1, 
n = 160000). 

The next, relatively long stage of evolution is char- 
acterized by restless moving and strong competition of 
RSWs. As the matter of fact, this stage is of the "winner 
takes all" (WTA) type. So, at the end of this stage, the 
single, fully regular RSW with Qt = —I occupies the 
whole active medium (Figure 1, n = 365000). 

But the core of the winner is not far from the active 
medium boundary. Drift of this RSW leads to collision 
of its core with boundary, and the core is absorbed by 
the latter (Figure 1, n = 391000). The residual non- 
spiral autowaves gradually run out the active medium 
(Figure 1, n — 391000) and excitation collapses fully at 
n — nc — 396982 (not shown at Figure 1). 

The described scenario of excitations collapse is very 
typical but not the single one. Some other ways lead- 
ing to collapse will be described in a forthcoming pub- 
lications. Now we will consider cases of non-collapsing 
systems. 



B. Dynamic Stabilization, Synchronization and 
Coexistence of RSWs 

In this Subsection, we describe a typical scenario of the 
2D-3L-2C excitable system evolution leading to dynamic 
stabilization and coexistence of RSWs. In this case the 
system reaches one of its periodic spatio-temporal attrac- 
tors — fully synchronized cyclic transitions (Li Lm 
Lii ^ ii ^ ...) for all AGs. Some typical stages of such 
the scenario of evolution are shown at Figure 2. 

The system evolves at first stages (0 < n < 10"*) to a 
mixed RSW-labyrinthic structure (Figure 2, n = 10000), 
where nucleation of RSWs has more pronounced form, 
than at Figure 1 (for the same n = 10000. 

All the subsequent stages of evolution at Figure 2 dif- 
fer more and more comparatively to Figure 1 . Instead of 
formation of a giant vortex, which pushes the labyrinthic 
structure to the boundaries of the active medium (Fig- 
ure 1), at Figure 2 we see multiple emergent RSWs 
[n = 30000). These RSWs have different sgn(QT), but 
RSWs with sgn((3T) = —1 are prevailed here (Figure 2, 
n = 30000). 

As the result of such direction of evolution, at the next 
stage only negative-charged, well-formed RSWs with 
Qt = — 1 occupy the whole active medium (Figure 2, 
n — 100000) — instead of the single giant multi-charged 
vortex with IQtI > 1 (Figure 1, n = 100000). 

Due to complex interactions between these RSWs as 
well as between RSWs and boundaries, the processes of 
revival of positive-charged RSWs take place (Figure 2, 
n = 165000). In contrast to Figure 1, where strong com- 
petition between RSWs leads to WTA dynamics, compe- 
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tition is weakened in the case under consideration (Fig- 
ure 2, n = 204000). 

Weakening of competition and low mobility of RSWs 
determine qualitatively another final (Figure 2, n = 
300000 and n — 987000) of evolution scenario (compara- 
tively to Figure 1). Namely, the dynamic stabilization of 
RSWs over the active medium is observed at n « 300000. 
Configuration of the system is precisely repeated (com- 
pare n = 300000 and n = 987000 at Figure 2). In other 
words, this is spatio-temporal limit cycle. 

Fully synchronized cyclic transitions for all AGs (Li 
Liu Lii Li ...) for this system (Figure 2) is in 
contrast to collapse of excitations for the previous sys- 
tem (Figure 1) despite of obvious transient ordering of 
the spatio-temporal dynamics in the course of events for 
both the cases (Figure 1 and Figure 2). Note that in 
both these cases we deal with spontaneous decreasing of 
effective freedom degrees quantity during the evolution 
of the system . 

Transient ordering of this kind, leading to collapse 
of excitations (Figure 1), may be interpreted as self- 
organized degradation (self-destruction) of the system. 
Final state of such the system, of course, is not self- 
organized, but it is still highly ordered. This kind of 
static order is known as freezing, because all AGs of the 
system occupies their ground levels (like usual particles 
of a physical system at absolute zero temperature) . 

In contrast, transient ordering, leading to dynamical 
stabilization of more or less complex spatio-temporal 
structures (Figure 2), may be interpreted as "normal" 
self-organization in autonomous dissipative system 
when most degrees of freedom are "slaved" by the small 
rest of ones. At this self-organized state, the system 
reaches the end of its evolution. But the system remains 
unfrozen, it is still far from equilibrium and possesses 
obvious order in the form of spatio-temporal limit cycle 
(Figure 2). All parts of this system (namely RSWs) are 
well defined and robust (as individual dynamical spatio- 
temporal autostructures). At the same time, all these 
parts and entirely synchronized at global level (as dy- 
namical components of the whole system demonstrating 
collective rhythmic motions). 

Nevertheless, both these scenarios of dissipative sys- 
tem evolution (self-destruction at Figure 1 and normal 
self-organization at Figure 2) exhibit more or less fluent 
transition from initial spatial chaos to static (Figure 1) 
or dynamically stable (Figure 2) order. 

A different scenario is exemplified in the next Subsec- 
tion. 



C. Transient Chaos Caused by Competition of 
RSWs and Labyrinths 

In this Subsection, we will describe another way to 
reach periodic state in a bounded discrete 2D-3L-2G ex- 
citable system. Formally, the final state for this new case 
is spatio-temporal limit cycle too (as in the case of Fig- 



ure 2). But the time of transient process is much longer 
here, and the longest part of the transient proceeds in 
the form of lethargic competition between various spatio- 
temporal structures. 

As it was already shown in the preceding Subsection 
(Figure 2), spatio-temporal structures may actively in- 
teract not only between themselves, but they interact 
in a nontrivial manner with boundaries of the excitable 
medium too. The simplest case when drift of RSWs leads 
to collisions of their cores with boundary, and the cores 
are absorbed by the latter (Figure 1, n = 391000), is 
not typical case for the system studied in this Subsec- 
tion. On the contrary, RSWs in such the system under 
appropriate conditions may be "refiected" from bound- 
aries in a strongly nonlinear manner^ . Moreover, due 
to nonlinear processes in the surface layer, the quantity 
of "reflected" RSWs may exceed the quantity of primary 
RSWs (phenomenon of multiplication of RSWs). A de- 
tailed description of nonlinear " reflections" of RSWs and 
some accompanied phenomena will be published in sep- 
arate paper. 

Here we describe a particular, but important case of 
highly bottlenecked transient process, caused mainly by 
almost everlasting competition of RSWs and labyrinths. 
At Figure 3, a fragment of evolution is shown at n > 10^ 
for the 2D-3L-2G system, which differs from the previous 
system (Figure 2) only by two parameters: = 44; h — 
46 for Figure 3 (instead of = 46; ft- = 51 for Figure 2). 
Starting pattern for the system at Figure 3 is precisely 
the same as at Figure 2 (so it is not shown at Figure 3). 

Previous two systems (Figures 1 and 2) reach their at- 
tractors already at n < 5 • 10^. The system under consid- 
eration (Figure 3) has very long transient time, because 
of lethargic competition between RSWs and labyrinthic 
structures under conditions of regeneration (nonlinear 
"reflections") and multiphcation of RSWs at boundaries. 
There are many fine phenomena accompanied these com- 
petition. E. g., the largest RSW with Qt = +1 resided 
in bottom left corner of the active medium (Figure 3, 
n — 1307000) moves in the North direction to the big 
labyrinth. Finally, the latter absorbs the core of the 
pointed RSW (Figure 3, n = 1322000). At the same 
time new small RSWs are generated in the bulk active 
medium and regenerated (by nonlinear "reflections") in 
the surface layer. 



^ Reflection of RSW from boundary, as well as reflection of any 
other moving dissipative structure (autostructure, autowave) is 
usually considered as forbidden. This is correct, commonly 
speaking, only in lowest (weakly nonlinear) approximations, 
when the properties of surface layer are assumed as close to prop- 
erties of bulk active medium. But real processes in surface layer 
may lead to strong nonlinear phenomena of, e. g., revival (re- 
generation) of RSWs or even to their duplication, triplication 
etc. Complex phenomena of such the kind may be considered 
as higher-order nonlinear interactions of RSWs with boundary, 
even if boundary conditions itself are "simple" (e. g. zero-flow 
conditions, as in the present work). 
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The dominant drift directions of RSWs in the system 
under consideration He at North^South and East^ West 
Unes (this circumstance is evident when one see a movie 
of the system evolution). Drift of labyrinths is more 
complicated. In the fragment shown at Figure 3, the 
big labyrinth is slowly moving approximately in North- 
East direction. A new labyrinth is formed simulta- 
neously. This new labyrinth pushes RSWs (Figure 3, 
n = 1322000) etc. 

The fragment of evolution shown at Figure 3 demon- 
strate only small part of transient phenomena observed 
during this system evolution. We will mention now such 
the typical stages of the system evolution {not shown? at 
Figure 3): 

• Nonlinear "reflection" of large RSW with \Qt — 1 
from boundary accompanied by changing of sgn((3T) and 
simultaneuos birth of several new small RSWs (satellites) 
nearby the core of large RSW (672000 < n < 687000). 

• Rotation of giant labyrinth around the domain oc- 
cupied by RSWs (2060000 < n < 2150000); erosion of 
this labyrinth due to RSWs activity up to almost full 
dominance of RSWs over active medium (2160000 < 
n < 2180000); revival of labyrinth structures and re- 
sumption of competition between RSWs and labyrinths 
(2190000 <n< 2210000). 

• Nonlinear "reflection" of large RSW with \Qt = 1 
from boundary with saving of sgn^Qx) (2309000 < n < 
2324000, upper part of active medium). Shining exam- 
ples of multiplication of RSWs at boundaries (the same 
interval 2309000 <n< 2324000, left, right and bottom 
boundaries). 

• Birth of self-organized pacemaker (a source of con- 
centric autowaves) due to nonlinear processes of RSW's 
core transformations in the surface layer; evolution and 
decay of the pacemaker (2965000 < n < 2980000, left 
boundary of the active medium). 

• Pushing of labyrinth by a group of RSWs (3100000 < 
n < 3115000) etc. 

From the formal point of view, all these stages of the 
system evolution are only partial transient episodes on 
the long way to the final self-organized state. We know 
that such the state is a regular attractor (spatio-temporal 
limit cycle or at least point-like attractor) because by def- 
inition our system is autonomous, it has bounded quan- 
tity of discrete elements, each element has bounded quan- 
tity of discrete states, and time is discrete too. But we 
do not get concrete view of the attractor until we reach 
it. Shortly, "transient is nothing, attractor is all". 

From an alternative point of view, each stage of the 
system evolution is a day in the life of this system, and 
the final of evolution is of minor interest, because no 
new events will come to pass if an attractor is already 
reached. Shortly, from this point of view, "transient is 



^ Additional figures to this Preprint may be requested from the 
author by e-mail (see the first page of the Preprint) 



all, attractor is nothing" . 

But this alternative is not the single one. Another 
alternative is as follows: "transient is all that we can see, 
because attractor is unachievable" . 

This last alternative (for a non-CA system with con- 
tinuous spectrum of states of its elementary components) 
was discussed in a seminal work of J. P. Crutchfield and 
K. Kaneko 6J, which was entitled: "Are Attractors Rel- 
evant to Turbulence?" . In , this question was slightly 
reformulated in the context of our studies of CA systems: 
"Are Attractors Relevant to Transient Spatio- Temporal 
Chaos?" (turbulence in bounded, fully discrete system 
is no more than metaphor). The answer is "Yes" if an 
attractor may be reached for a reasonable time Tattr (m 
fully discrete and bounded system Tattr is always limited 
by quantity of all possible states of the system). But 
the answer is "No" if Tattr exceeds any possible duration 
of an experiment. In this case a system with transient 
spatio-temporal chaos cannot be distinguished from true 
chaotic system without additional testing. 

As a matter of fact, there are some intermediate classes 
of phenomena "at the edge between order and chaos" 
which may appear in bounded discrete deterministic 
system with large phase space. Self-organization sce- 
nario which includes super-slow, bottlenecked, chaotic- 
like stages is a signature of dominance of such an interme- 
diate class of system dynamics in numerical experiments. 
Really, having limited time and computer capacity, one 
cannot reach final self-organized state for a system with 
huge dimension of phase space. Such CA (or another 
discrete mapping with lethargic evolution) does not per- 
mit direct forecasting of the system future without direct 
computation. So the computed part of transient process 
is the single source of available information of our fully 
deterministic but partially determined system. 

Generally, the ideas of interconnection between chao- 
tization, turbulence, unpredictability (at one side) and 
ordering, self-organization, long-time forecasting (at the 
opposite side) may be very fruitful at least for investiga- 
tion of complex deterministic systems. 



V. CONCLUSIONS 

In this work, we fulfill computer modeling of spatio- 
temporal dynamics in large dicrete systems of three- 
level excitable ACs interacting by short-range 2C diffu- 
sion. Computer experiments with this 2D-3L-2C sys- 
tem were carried out using the cross-platform soft- 
ware package "Three-Level Model of Excitable System" 
(TLM) © 2006 S. D. Makovetskiy 

The most typical scenario of evolution is as follows. A 
robust RSW or more complicated but strictly periodic in 
time structure is formed being a cyclic attractor at sev- 
eral initial conditions — see Figure 2. This is an analog 
of limit cycle, i. e. regular attractor known in lumped 
dynamical systems. Alternatively, collapse of excitations, 
i. e. full freezing of the excitable system takes place — 
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sec Figure 1. This is an analog of point-like attractor in 
lumped dynamical systems. 

Long-time evolution (10*' — 10^ iterations) of a 2D- 
3L-2C system with slightly changed parameters demon- 
strates some unusual phenomena including highly bottle- 
necked collective relaxation of excitations. Part of these 
phenomena are caused by mechanism of nonlinear regen- 
eration of RSWs in the surface layer of an active medium. 
Instead of well known scenario of RSW collapse, which 
takes place after collision of RSWs core with absorb- 
ing boundary, we observed complicated transformations 
of the RSWs cores leading to nonlinear "reflections" of 
the RSWs from the normally absorbing boundaries or 
even to birth of new RSWs in the surface layer. Lethar- 
gic transient processes observed in our computer experi- 
ments are partially caused by competitions between per- 
manently reviving RSWs and labyrinthic spatio-temporal 
structures. To our knowledge, phenomena of such the 
nonlinear "reflections" of RSWs and resulting die hard 
vorticity in excitable media with absorbing boundaries 



were unknown earlier. 

The author is grateful to D. N. Makovetskii (Institute 
of Radio-Physics and Electronics, National Academy of 
Sciences of Ukraine) for stimulating discussions on vari- 
ous aspects of this work. 



APPENDIX A: LIST OF ABBREVIATIONS 

IC, 2C, ... One-Channel, Two-Channel, ... 

ID, 2D, ... — One-Dimensional, Two-Dimensional, ... 
2L, 3L, ... — Two-Level, Three-Level, ... 

AC — Active Center 

CA — Cellular Automaton 

CP — Control Parameters 

GBC — Geometry and Boundary Conditions 

RSW — Rotating Spiral Wave 

TLM — Three-Level Model (of excitable system) 

ZM — Zykov-Mikhailov 
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FIGURE CAPTIONS 

to the paper of S. D. Makovetskiy 

"Numerical Modeling of Coexistence, Competition and Collapse 

of Rotating Spiral Waves in Threc-Lcvcl Excitable Media 
with Discrete Active Centers and Absorbing Boundaries" 

(figures see as separate PNG-files) 

Figure 1: Birth, evolution and collapse of RSWs with the same (by magnitude and sign) effective 
topological charges. Dimensions of the active medium are Mx = My = 300. Black, gray and white pixels denote 
ACs in excited (im), refractory (in), and ground (Li) states respectively. Starting pattern is shown at n = 0. The 
relaxation times of ACs = 50. The set of CPs is as follows: g = 1; h = 50; / = 3. 

Figure 2: Birth, evolution, dynamic stabilization and coexistence of RSWs. Dimensions of the active 
medium and colors of pixels are the same as at Figure 1. Starting pattern (n = 0) is not the same as at Figure 1, 
but statistical properties of both patterns are almost identical. The relaxation times of ACs and the set of control 
parameters are as follows: Tg = 60; Tr = 4:6; g = 1; h = 51; / = 2. 

Figure 3: Typical stage of transient spatio-temporal chaos (competition of RSWs and labyrinthic 
structures) at n ^ max(Te, Tr). Dimensions of the active medium and colors of pixels are the same as at Figures 1 
and 2. Starting pattern is precisely the same as at Figure 2 (and not shown here). The relaxation times of ACs and 
the set of CPs are as follows: Tg = 60; Tr = 44; 51 = 1; /i = 46; / = 2. 
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